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The shear and bulk viscosities of gluon plasma are calculated by accumulating a large amount of 
data for the Matsubara Green function (Gp(t m )) on isotropic 24 3 x 8 and 16 3 x 8 lattices. In the 
case of Iwasaki's improved action, the calculations of Gp (?„,) are carried out on roughly 6 million 
configurations, while for the standard action the calculations are done on more than 16 million 
configurations. The shear viscosities increase roughly with T 3 , and rj/s ratios are close to the 
KSS lower bound in the region where 1 < T /T c < 25. Using these data the bulk viscosities are 
also determined in the region where T/T c < 2. They are roughly one order of magnitude smaller 
than the shear viscosities. Our next target is to determine the transport coefficients more precisely 
by a maximum-entropy method. For this purpose the most effective method may be to adopt an 
anisotropic lattice. In this report, we study the possible systematic error due to the deformation 
of the anisotropic lattice at short distances. Near the critical temperature, it is found that the 
standard action suffers from a large deformation on the anisotropic lattice at short distances, 
while the deformation is slight for Iwasaki's improved action. To reduce the fluctuation of the 
Matsubara Green function, the improvement of the energy momentum tensor operator by using 
clover-type loops is promising. We are also attempting to apply the multi-level algorithm to 
reduce fluctuation. 
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1. Introduction and formalism 

A new state of matter has been reported in RHIC experiments. From an estimation of the 
temperature, it is expected that quark gluon plasma (QGP) is realized there. A phenomenological 
study of elliptic flow revealed that a new state of matter is well explained by the fluid with very 
small viscosity. Jet quenching data indicate that it is strongly interacting. 

If QGP is realized as a fluid, its viscosities should be smaller than those of usual fluids such 
as water and oil, because they are proportional to g~ 2 , and the coupling constant g is larger than 
1/137 even above the transition temperature. The lower bound for the ratio of shear viscosity(rj) 
to entropy(s) has recently been conjectured by Kovtun et al. [JT|] . However, values close to transi- 
tion temperature should be calculated from QCD by fully taking into account the nonperturbative 
effects. It is urgently necessary to carry out detailed phenomenological studies of the new state of 
matter using these viscosities. 

The transport coefficients are calculated in the framework of linear response theory: 

r] = - f d 3 x' f She E{t ^ f" dt' < T l2 {x,t)T l2 (x',t') > ret (1.1) 

J J— CO ,7 CO 

^77 + £ = -J d 3 xf J dt ie £ ^ p dt' < T n (x,t)T n ',t') > ret (1.2) 

where tj is shear viscosity, and £ is bulk viscosity. < T^ v T pa > ret is the retarded Green function of 
the energy momentum tensor. In pure gauge theory it is given by 

Tfiv = 2Tr[Ff ia F va — — SfiyFpaFpa], (1-3) 

and the field strength tensor F^ v is defined by the plaquette operator on a lattice: 
U^ v (x) =exp(ia 2 gF^ v (x)). 

The transport coefficients are also expressed by the slope of the spectral function p(co) at 
ft) = of the corresponding retarded Green function. The shear viscosity v\ is written as 

p(co) 

77 = 71 lim — — - . (1.4) 
m^O ft) 

On a lattice, we determine p(ft)) by the Matsubara Green function Gp instead of the retarded 
Green functions themselves, because both Green functions have the same spectral function, 

„ / s 1 v imt r P( w ) , f°° cosh(co(t m -B/2)) . „ cs 

G B (t m ) = --Y e^ m '''" / HK ' dco= , n ir\ \ p(co)dco, (1.5) 

p /3Y J-ooicOn-co Jo sinh(cop/2) 

where co n = 27in/f5 for the Matsubara Green function. If /ft),, is replaced by po + is, the retarded 
Green function is obtained. 

However, there are still difficulties in the determination of p (ft)). One is that Gp (t m ) is discrete, 
while p(ft)) is continuous. Therefore, fine resolution in the temperature direction (simulation on 
large Nt lattice) is necessary for its accurate determination. The other difficulty is that Gp(t) 
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is noisy; thus its determination requires much CPU time. Hence, we start with the smaller Nt 
lattice, assuming a plausible form for the spectral function p (ft)) that fits Gp (t m ) well. The simplest 
nontrivial form is|§], 

Ay y 

p(ftl)=-(- ^ ^ y) (1.6) 

KV ' % y (m - oo) 2 + y 2 (m + ft>) 2 + y 2/ 
This form is derived from a perturbative calculation in <p 4 theory [||]. 



2. Numerical results from isotropic lattice 

Because p(ft)) given by Eq.([OJ) has three free parameters, Gp(t m ) should be calculated on 
Nt>& latticels. We carry out simulations on 24 3 x 8 and 16 3 x 8 lattices, using Iwasaki's improved 
action and the standard action, and the temperature range is 1.4 < T/T c < 25. We are attempting 
to overcome the huge fluctuations by a large number of measurements. In the case of the improved 
action, Ga(t m ) is determined by roughly 6 million measurements, while for the standard action 
there are more than 16 million measurements. The fit of Gp(t m ) is made by SALS, and errors 
are estimated by the jackknife method. Then the viscosities are obtained by the formula T]a 3 = 
4Ayra/(y 2 + m 2 ) 2 . 

To obtain the viscosities in physical units, we need the lattice spacing a(g). For the improved 
action, a has been determined for 2.2 < j3 < 3.8 by the Tsukuba group[Q], and for the standard 
action, a has been determined for 5.58 < j3 < 6.5 by Edward et al.[Q]. Outside these regions, we 
assume a two-loop asymptotic scaling relation. 

The results for shear viscosity tj in physical units are shown in Fig.[jJ Because the t]xa 3 
has weak T dependence, the shear viscosity increase roughly with T 3 throughout the temperature 
region. We find little difference between the results from 24 3 x 8 and 16 3 x 8 lattices. Thus, the size 
effect may not be large for the lattices of these sizes. However, more accurate data are necessary to 
determine quantitative size dependences. 

We have also shown the perturbative results beyond the leading log approximation^] in Fig.[j], 
where the scale factor p in the running coupling constant is a free parameter. The agreement 
improves when p becomes smaller, but in this case the breakdown of the perturbative calculation 
starts at a higher temperature. 

Let us proceed to the r]/s ratio, recently studied by Kovtun et al.[jlj]. The results are shown 
in Fig.||[ The lattice data on entropies are given by Ref. [Q 0]. The perturbative results are also 
shown, where we use the entropy calculated by the hard thermal loop approximation given in 
Ref.[|8|]. Because both Tja 3 and so? have weak T dependence in the T > 1.5!T C region, the tj/s ratio 
also has weak T dependence and the lattice results are close to the KSS bound for T /T c < 25. 

In Fig.|3| we show the ratio in RHIC temperature regions together with the bounds by Meyer 
[Q], who employed a multi-level algorithm and got bounds without using an ansatz. 

For the improved action, signals for the bulk viscosities begin to overcome the errors, when 
there are about 6 million measurements. The results are shown in Fig.|[ The values of bulk 
viscosity still have rather large errors, but at T close to T c , their values are determined, which 
are roughly one order of magnitude smaller than the shear viscosities. Their T-dependence is an 
interesting problem but that requires more measurements. Recently it has been claimed that the 
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Figure 1: Shear viscosity in physical units from lattice and perturbative calculations. The circles are results 
from the improved action and squares are those from the standard action. The perturbative results beyond 
leading log approximations Q| are shown by lines. 




Figure 2: rj /* obtained by lattice simulations (circles and squares) and perturbative calculations (lines). The 
KSS bound[|IJ is also shown. 
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Figure 3: Ratio of 77/s in RHIC regions. 



bulk viscosity is large near the critical temperature and decreases rapidly with TjlOj]. Our results 
do not contradict this. In the case of the standard action, the bulk viscosities still have large errors 
that they cannot be determined. 

It will be interesting to carry out phenomenological studies on RHIC data taking into account 
these viscosities in the fluid model. 



2.1 Discussions 

• The renormalization factor Z of the energy momentum tensor is discussed by Meyer^]: 
Z = 1 — g 2 /2{c a — c T ). If the parametrization of Z factor given by Ref.[||] is used, the viscosities 
calculated by the standard action decrease by about 30%. Z factor can also be written as follows: 



Z = dy/dE,\ 14], where t, is the renormalized anisotropy and 7 is the bare anisotropy. In the case 
of Iwasaki's improved action, £ ~ 7 over a wide range of /3 and E, ; therefore the Z factor is close 
to 1. If the Z factor is taken into account, the difference between 77 obtained from improved action 
and standard action decreases. 

• We have attempted to fit G^ 2 by other parametrizations of p (ft)) than that given in Eq.[0]. If 



we apply the formula for p proposed in Ref.[ll], the fit is not satisfactory and p does not satisfy 



the constraint cop(co) > 0[12]. If we truncate the Taylor expansion of p(co) after the lowest 3 
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Figure 4: Lattice results of £/s ratio from the improved action 

terms, the fit is also not satisfactory and the coefficient of co becomes negative, which also cannot 
be accepted as a spectral function. In the three-parameter functions for p, we have not found a 
p((o) that fits Gg 2 (? m ) well except for EqjL6|. 

3. Next project: toward the high-precision calculation of r\ 

3.1 Simulation on anisotropic lattice 

To determine p(ffl) without relying on any assumption, the maximum entropy method(MEM) 
seems promising. To get a reliable result, accurate Gp for a Nj ~ 30 lattice is necessary. For this 
purpose, the best method may be to adopt an anisotropic lattice. Before carrying out the simulation 
on large anisotropic lattices, we start with a preliminary study on the possible systematic error. 

The fundamental properties of anisotropic lattices have already been studied[|l3|, 14, 15]. The 
anisotropy £ = a a /a T is controlled by the bare anisotropy y in the action. The y dependence of % 
is expressed in terms of the ratio E, / /JO], |l4[ |L5] ]. The £/y ratio is determined by its asymptotic 
plateau for r > 3, where r is the distance of the lattice in the unit of space direction a a . At short 
distances, the t, //ratio is not equal to its asymptotic value. Thus Gn(t m ) for small t m /E, will suffer 
from a systematic error due to deformation. 

To study the effects of the deformation, we compare G 1 ^ for isotropic and E, = 2 anisotropic 
lattices at the same T/T c (a a ). For the standard action, G 1 ^ is shown in Fig.|[ Large discrepancies 
are observed near T c , due to the deformation at short distances (t m / § < 2) and to the difference 
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Figure 5: G» 2 (f m ) on isotropic (iVr = 4) and <^ = 2 anisotropic lattice (Nt = 8) 



between the Z factor in the lattices. An increase in the deformation is anticipated as the E, /y ratio 
increases toward T c . On the other hand in the case of improved action, the difference in G l p(t m ) for 
the lattices is small. In this action, because the ratio t, /y is close to unity over a wide range of jS 
and t, , the deformation is weak and the Z factor is close to unity. 

3.2 Improvement of the energy momentum operator 

To reduce the fluctuation of Gp , we have attempted to use clover-type operator for the defini- 
tion of the energy momentum tensor. It was found that the noise is strongly suppressed, and the 
normalization and the ?-dependence of G^ v (t) are not changed. We are also testing the effective- 
ness of the multi-level algorithm^] for reducing the fluctuations[17]. 



4. Conclusion 

The shear and bulk viscosities are calculated on isotropic 24 3 x 8 and 16 3 x 8 lattices. The 
T-dependence of the t]/s ratio is weak and its value is close to KSS bound throughout the region 
T /T c < 25. The bulk viscosities are obtained in the T /T c < 2 region, and their values are one order 
of magnitude smaller than those of shear viscosities. 

The systematic error due to the deformation of lattice spacing at a short distance is studied. 
For the standard action, the deformation becomes strong near T c . In the t m /E, < 3 region, the effects 
of deformation should be carefully controlled. 

For the accurate determination of the spectral function p by MEM, improvements of the energy 
momentum tensor operator using clover-type loops and the multi-level algorithm are promising. 
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